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ABSTRACT 

Motivated by the recent detection of stochastically excited modes in the massive star V1449 Aql (Belkacem et al., 2009b), already 
known to be a f} Cephei, we theoretically investigate the driving by turbulent convection. By using a full non-adiabatic computa- 
tion of the damping rates, together with a computation of the energy injection rates, we provide an estimate of the amplitudes of 
modes excited by both the convective region induced by the iron opacity bump and the convective core. Despite uncertainties in the 
dynamical properties of such convective regions, we demonstrate that both are able to efficiently excite p modes above the CoRoT 
observational threshold and the solar amplitudes. In addition, we emphasise the potential asteroseismic diagnostics provided by each 
convective region, which we hope will help to identify the one responsible for solar-like oscillations, and to give constraints on this 
convective zone. A forthcoming work will be dedicated to an extended investigation of the likelihood of solar-like oscillations across 
the Hertzsprung-Russell diagram. 



1. Introduction 

/? Cephei-ty pe stars are known to p ulsate with large amplitude 
oscillations ( bautschv & Saioll 1995L[l996h . The Ar-mechanism is 
at the origin of those large amplitude oscilla tions and is related t o 
the iron opacity bump at log 7/ m 5.3 (e.g. iPamvatnvkhl 1 1 999b . 
Those modes are linearly unstable. In contrast, low-amplitude 
modes have been observed in the Sun for years and are thought 
to be intrinsically stable but stochastically excited by turbulent 
convection. Their excitation is attributed to turbulent convection 
and takes place in the uppermost part of the Sun, which is a place 
of vigourous and t urbulent motions. Since the pioneering work 
of iLighthillld 19521) . we know that a turbulent medium generates 
incoherent acoustic pressure fluctuations (also called acoustic 
"noise"). In the last decade, solar-like oscillations have been de- 
tected in numerous F-G type main sequence and red giant stars, 
in different evolutionary stages and with different m etallicities 
(see the recent review bv lBedding & Kieldsenll2007l) . 

As in the Sun, massive stars have convective regions, namely 
the convective core, and two regions associated with the helium 
and iron opacity bumps. The inner convective region (i.e. the 
convective core) has an important impact on the internal struc- 
ture and on the subsequent evolution of the star through, for 
instance, th e extent of the region where c hemical elements are 
mixed (e.g. iKippenhahn & Weigerll 1 1990h . In contrast, the ex- 
ternal convective zones associated with opacity bumps are gen- 
erally thought to be unimportant since they transport only a very 
small fraction of the energy flux, which is domi nated by radia- 
tion. It is only since the revision of the opacities dlglesias et al.L 
which led to an enhancement of the iron opacity bump 
and could thus account for the instability in fi Cephei stars, that 
a convective zone related to the iron opacity bump has been 
thought to be present in massive stars. Much attention is now 



drawn to this convective region associated with the iron opac- 
ity bumpQ, because it is of interest for the understanding of 
surface effects such a s microturbulence or wind clumping (e.g. 
ICantiello et all l2009t) . However, both the convective core and 
the iron convection zone are poorly understood. In particular, the 
typical length-scale and time-scale of turbulent motions associ- 
ated with these regions are only inferred through mixing-length 
arguments, as well as dynamical properties such as the convec- 
tive velocities. 

Solar-like oscillations (linearly damped and stochastically 
excited modes) are closely related to those dynamical prop- 
erties, so their detection in such massive stars is a promising 
way to infer them. Until recently, stochastically excited modes 
have only bee n predicted and detec t ed in s olar-like and red gi- 
ant stars (see iBedding & Kieldsenl 120071: iMichel et all 120081: 
Ide Ridder et al.L 120091 for details). The data gathered by the 
CoRoT mission allowed us to report the first d etection of solar 
like os cillations in a massive star, V1449 Aql (iBelkacem et all 
l2009bl) . Those oscillations exhibit mode amplitudes higher 
than solar, hence the issue of their excitation is of inter- 
est. Investigation of unstable modes already gives constraints 
on the internal structure of massive star s (e.g. [Aertsetal , 



| 2003 | :lThoul et all 120031: Ibupretet dl l2004t iPamvatnvkh et al 
3041) 



2004 ) as well as on the underlying excitation mechanisms (e.g . 
Moskalik & Dziembowskil, |1992[ Dziembowski & Pamiatnvkh , 
1993c iDziembowski et al l ll993tlDziembowski & Pamvatn vkh, 



2008). The discovery of solar-like oscillations is thus an oppor- 



tunity to go a step further i n the understanding of those stars 
which still challenge theory dDziembowskil [2007b . 

Our objective in this work is to assess the excitation of 
solar-like oscillations in massive main-sequence stars, such as 
P Cephei type stars, by turbulent convection. We aim at deter- 
mining which convective region is able to excite modes above 
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1 In the following, this convective region will be referred to as iron 
convective region. 
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the observational threshold, as well as emphasising the possible 
asteroseismic diagnostic that would permit us to identify the ex- 
citation region and to obtain constraints on this region. As the 
amplitude of stochastically excited modes is a balance between 
driving and damping, both must be modelled. In this work we 
will use the formalism proposed bv | S a madi & Goupill(l200ll) and 
extended bv lBelkacem et all (2006a, b. 2008). Damping rates are 
comp uted using the MAD code dDuprettl2002l:lGrigahcene et all 
120051) . which is based on a full non-radial non-adiabatic treat- 
ment. This first paper of the series focuses on a bench- 
mark model. An extended investigation of the mode amplitude 
through the Hertzsprung-Russell diagram will be the object of 
the second paper of this series. 

The paper is organised as follows: Sect. 2 introduces the 
general formalism that is used to compute mode amplitudes. In 
Sect. 3, we focus on mode damping rates and we emphasise their 
frequency dependence. Excitation by the convective region asso- 
ciated with the iron convective region is also studied and a pos- 
sible diagnostic is inferred. Section 4 presents the computation 
of mode amplitude excited by the convective core. Conclusions 
are formulated in Sect. 5. 



2. Modelling mode amplitudes 

2. 1 . Equilibrium stellar structure; a 10 M Q model as a 
benchmark 

The objective is to explore the potential of convective regions 
to drive oscillations. As a benchmark we investigate the driv- 
ing of radial acoustic modes by turbulent convection zones in 
a 10 M main-seque nce star, which is typical of j3 Cephei stars 
dPamvatnvkhl U999). A more systematic exploration of the HR 
diagram is dedicated to the second paper of this series. 



2.1 .1 . Physical input 

The stellar structure models used in this work are obtain ed with 
the stellar evolution code CLES (IScuflaire et all I2008I). In the 
interi or model, we used the OPAL opacities ( Iglesias & RogersL 
1996) extended to low tem peratures with the opacities of 



Alexander & Ferguson! d 19941) and the C EFF equation of state 
dChristensen-Dalsgaard & Dapper! fl 9921) . The metallicity is as- 
sumed to be Z = 0.02 with the metal mixture as derived by 
iGrevesse & Noels! (119931) . Convection is included according to 
a Bohm-Vitense mixing-length (MET) formalism and no over- 
shoot is allowed. The effect of the mixing-length parameter is 
discussed in Sects. [331 and 14. 3| Note, however, that the default 
value of the mixing-length is 1 .8 (the Solar value) and in the con- 
vective region associated to the iron opacity bump, the mixing- 
length is of the same order that the extent of the convective re- 
gion, 1.5% of the star radius. 

The location of the benchmark model on the Hertzsprung- 
Russell diagram is plotted in Fig.Q] The model has an effective 
temperature of log r e ff = 4.35, and a surface gravity of logg = 
3.80. 

2.1.2. Convective zones 

In such massive stars, one can distinguish three convective re- 
gions, namely; 

- the convective zone associated with the helium II opacity 
bump. This region is located at T ~ 40 000 K (Pig. [2j, and is 
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Fig. 1. Hertzsprung-Russell diagram for a ten solar mass model 
computed with the CLES code (see Sect. 12. U . The diamond cor- 
responds to the benchmark model as described in Sect. l2.1.Tl the 
triangle symbols corresponds to another model used for com- 
parison in Sect. 13.2.21 Star symbols correspond to the different 
evolutionary stages displayed in Fig. [6] (bottom). Eventually, the 
square and triangle symbols are used to identified the models 
presented in Fig.|6](top). 



very weak in the sense that it is inefficient in transporting en- 
ergy with a negligible ratio of the convective heat flux to the 
radiative one. In addition, such a region is located near the 
stellar surface, where the density is low, and velocities are in 
the 10 to 20 m.s~' range (low heat capacity is the cause of ef- 
fective smoothing of the temperature contrast and reduction 
of the buoyant acceleration). Hence, the kinetic energy that 
could feed modes is small. 

- the convective zone associated with the iron opacity bump, 
located at T « 200 00Qg (Fig. It has been shown 
dCantiello et all l2009t) that the occurrence of convection in 
the iron opacity peak is strongly dependent on the luminos- 
ity and the metallicity of the star. At a given luminosity, a 
smaller metallicity implies a lower opacity, which means that 
a larger fraction of the total flux can be transported by radia- 
tion. On the other hand, increasing the luminosity at a given 
Z leads to a more important contribution of convection in the 
energy transport. For a given metallicity, there is a luminos- 
ity threshold (of the order of 1O 3 2 L f° r Z = 0.02) above 
which a convective zone is associated with the iron peak. 
As for the helium convective region, the transport of energy 
by convection is very low. However, this region is located 
deeper than the helium one, where density is higher, with 
convective velocities up to 5.5 km.s~'. The kinetic energy 
flux is still negligible compared to the total flux but the ratio 
of the kinetic energy flux of the iron convective zone to the 
helium one is * 8 x 10 7 . 

- the convective core. Convection is fully adiabatic (fig.|2]i, and 
a significant part of the total flux is transported by convec- 
tion. The kinetic energy is four order of magnitudes higher 
than for the iron convective region, due to high densities, 
with velocities around 300 m.s~'. 

We conclude that the convective zone associated with the helium 
opacity bump is unlikely to be able to efficiently excite solar-like 
oscillations. Indeed, it has been verified numerically that the re- 
sulting theoretical mode amplitudes are found well below the 



Belkacem et al.: Solar-like oscillations in massive main-sequence stars 



3 



> 
> 



0. 1 5 


...... .. .. 1 


... i . .. .. , 


0.10 




Helium C.C. ; 


0.05 




Iron C.C. / 1 J 


-0.00 


! Convective core 




-0.05 






-0.10 






-0.15 


.... i .... i .... i 





7.5 



7.0 



6.5 



6.0 
log T 



5.5 



5.0 



4.5 



Fig. 2. The difference between the temperature gradient and the 
adiabatic one (V - V a d) versus the temperature. In the convec- 
tive core, the difference V - vanishes since convection is 
efficient and transport almost all the energy flux. In contrast, in 
the outer convective regions, we see significant gradients which 
means that convective transport is not efficient there. 



observational detection threshold of CoRoT (around one ppm). 
In contrast, the convective zones associated with the iron opac- 
ity bump as well as with the convective core are more promis- 
ing. The available kinetic energy flux is important. However, to 
be efficient, mode excitation must fulfil other criteria, which are 
discussed in Sects. |3~2l and|4l 



adiabatic pulsation code MAD dDupretl 120021) . This code in- 
cl udes a time-dependent c onvection (TDC) treatment described 
in lGrigahcene et ail d2005l) . which allows us to take into account 
the role played by the variations of the convective flux, the turbu- 
lent pressure, and the dissipation rate of turbulent kinetic energy. 

The damping of acoustic modes is found to be dominated by 
the perturbation of the radiative flux. Note that, although crucial 
for the excitation, the upper convective regions do not affect the 
damping rates for stars we are considering. Indeed, convection is 
very inefficient for the transport of energy (compared to the ra- 
diative transport) and the feedback of the convective background 
on the pulsation remains small. We have numerically verified 
that the interaction between convection and pulsation does not 
affect damping rates as well as the effect of turbulent pressure 
and the dissipation rate of turbulent kinetic energy. 

Hence, the dominant contribution to the damping rate can be 
written 



1 



1 = 



2 uqI 



M 



Im 



^ T SS 
Po 



dm 



(3) 



where 6S is the perturbation of entropy, dp the perturbation of 
density, 7q the mean temperature, a>o the mode frequency, po the 
mean density, m the local mass, M the total mass, and the star 
stands for the complex conjugate. 

Keeping only the radial contribution of the radiative flux in 
the energy equation because of its dominant contribution, and 
neglecting the production of nuclear energy (<5e = 0), one gets 



icrTodS 



ASL r 

dm 



(4) 



2.2. Computation of mode amplitudes 

We compute the me an-squared surface velocity for each radial 
mode as described in lBelkacem et al.l d2006b, 2009a|), i.e. 



v 2 M = 



P 



2t]M 



(1) 



where Ai = I/\^ r (h)\ 2 denotes the mode mass, / the mode iner- 
tia, 77 the damping rate, P the energy injection rate, and h is the 
height in the stellar atmosphere. In this section, we consider the 
level of the photosphere h — R with R being the radius at the 
photosphere. 

The mode amplitude in terms o f intensity is then ded uced 
at the photosp here, as proposed by iDziemblowskil i 19771) and 
lPesnelll(ll99C 



i.e. 



where cr - <±>q + it], and L r is the radiative flux. In addition, in 
the diffusion approximation, one can write 



5L 



5« | 1 dST | 2 fr | 3 6T 6p 8SA 
k (dr/dr) dr r T p dr J 



Finally, in the quasi-adiabatic approximation (which is valid in 
the deep layers) one has 



5k dp 
— SS K ps — with 
K p 



K ps = (k p + (T 3 - 1) k t ) IT\ (6) 



and 



dlnK 
d\np 



<91n/e 
dlnp 



,kt 



dlnT, 



(7) 



— = 1 STeS +2 — 
L T eS r 



(2) 



where 6L is the bolometric mode intensity fluctuation, and 6T e tf 
the mode effective temperature fluctuation. The relation between 
the variation of effective tem perature and the mode displacement 
is obtained from MAD (see iDuprell l2002t iDupret et all I2002L 
for details). Eventually, as the mode amplitude is a balance be- 
tween the driving (P) and the damping (77), both are to be deter- 
mined and computed. 

2.2.1. Computation of the damping rates 

As shown by Eq. (Q]), the computation of mode amplitudes re- 
quires the knowledge of both energy injection (P) and damp- 
ing (77) rates. The latter has been computed using the full non- 



2.2.2. Energy injection rates (P) 

The formalism we use to compute energy injection rates of radial 
m odes was developed bv[Samadi & Goupill(l2001l) and extended 
bv lBelkacem et all d2006allbl) . It takes two sources into account, 
the Reynolds stress tensor and the advection of the turbulent fluc- 
tuations of entropy by the turbulent motions (the "entropy source 
term"). The entropy contribution is found to be around ten per- 
cent of the Reynolds one. 

The energy injection rate, P, then mainly arises from the 
Reynolds stresses an d can be written as (see Eq. (21) of 
iBelkacem et~aTll2008l) 



P = 



16 7T 3 

1521 



dm po 



dr 



S R 



(8) 
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with 

X+°° ^k r+°° 
^ E \ k ) J dt0Xk(u + f>d)Xk(o) (9) 

where Sr is the source function, E(k) the spatial kinetic energy 
spectrum, Xk the eddy-time correlation function, and k the wave- 
number. 

The rate (P) at which energy is injected into a mode is 
then computed according to Eq. ([8]). Eigenfrequencies and eigen- 
functions are computed using the non-adiabatic pulsation code 
MAD. In addition to the eigenfunctions and the density stratifi- 
cation, Eq. ([8]l involves both the convective velocity and the tur- 
bulent kinetic energy spectrum. The typical convective length- 
scales are poorly known for main sequence massive stars. Hence, 
the classical mixing-length theory is used to find the injection 
length-scale (i.e. the scale at which the turbulent kinetic en- 
ergy &£ectrumjs_jn^imu^ and a parameter /3 is introduced 
(see lSamadi & Goupilll200il for details) such that the associated 
wave-number is ko = 2n//3 A, where A is the mixing-length. One 
has also to specify the way that turbulent eddies are temporally- 
correlated by defining an eddy-time correlation function. A 
Lorentzian function had successfully been used ; in the solar case 
dSamadi et al.L l2003bl iBelkacem et al.1 l2006bl). as we l l as fo r a 
Cen A, and HD49933 jSamadi et all 120081; ISamadl 120091) . it 
indeed reproduces the observational data. Consequently, if not 
explicitly mentionned, such modelling will be used in the fol- 
lowing. 

3. Driving by the iron convective region 

3.1. Damping rates 

3.1 .1 . Efficiency of the damping 

In Fig. [3] the fundamental mode is found linearly unstable with 
respec t to the /c-mechani sm associated with the iron opacity 
bump (iPamvatnvkhl 1 1999b . Higher radial-order modes are found 
linearly stable. The damping rates range between one and three 
/iHz for v > 150 /iHz, while they nearly vanish around the unsta- 
ble mode. They are found to be dominated by the first two terms 
ofEq.©. 

Conc erning the first contr ibution (5k/k), as described for in- 
stance by IPamvatnvkhl (1 1 999b . the region where k p + {Y?, - Y)kt 
increases outwards tends to drive the mode while the region 
where k p + (T3 - 1)kt decreases outward damps the mode. This 
contribution then mainly determines the damping rates of low- 
order modes while the second term of Eq. (O begins to be im- 
portant for higher order modes. 

3.1.2. Periodic variations 

For high frequency modes, the damping in the range log T e 
[5.2; 4.9] dominates over the driving region (log T e [5.5; 5.3]). 
This can be seen in Fig. |4] in which the work integral is plotted, 
for three modes which correspond to the first and the second 
maxima and the minimum in between of the damping rate (see 
Fig.0. 

The oscillation of the damping rates for high-frequency 
modes (v > 175 /iHz) is related to the location of the node of 
SP/P with respect to the damping region. By considering the 
pressure fluctuations (SP/P) of three modes, namely n = 7, 
n — 9 and n = 12, Figure |4] shows that the damping is mini- 
mum when there is a node in the damping region while maxi- 
mum when the nodes are located at the edges of those layers. 



3.0 p 




v [/-iHz] 

Fig. 3. Damping rates (77) of radial modes as a function of fre- 
quency, computed as described in Sect. 12.2.11 using the stellar 
model introduced in Sect. 12.11 The fundamental radial mode is 
found unstable (77 < 0) while the overtones are found stable 
(V > 0). 



The second maximum arises when two nodes of SP/P straddle 
the damping region. 

This is easily explained by the set of equations Eq. ([3])- 
Eq. since a node in pressure fluctuations lowers the opacity 
fluctuations that appear to be the main damping term in Eq. (0. 
Note also that the second derivative of ST I T is negative when 
ST /T is positive, which strengthens the radiative losses. 

It then creates the periodic variations observed in Fig. [3] for 
v > 175 //Hz. Hence, the detection of such an oscillation of mode 
line-width (damping) would be a signature of the iron opacity 
bump and also an opportunity to get constraints on modes pres- 
sure perturbation in this region. 

3.2. Energy injection rates 

The poor knowledge of the dynamical properties of convective 
regions associated with the iron opacity bump leads us to use a 
simple description based on the MET as explained in Sect. 12.2.21 
In terms of energy injection rates, two major features are to 
be determined: the convective velocity and the injection length- 
scale. As for the superficial convective layers of the Sun, the iron 
convection region is inefficient since it transports a negligible 
part of the energy flux through convection. The Mach number 
(Ai ~ 0.1) is similar to that found in the upper convective layers 
of the Sun in which M ~ 0.3. To what extent one can use the 
values of the parameters obtained for the Sun is however diffi- 
cult to assess and only numerical simulations will provide a firm 
answer. By default, for the iron convective region, we will use 
the values of the convective velocity provided by MLT using the 
assumed solar value (a = 1.8). In addition, the injection length 
scale is deduced from the numerical simul ations of the upper par t 
of the solar convection zone as derived bv lSamadi et al.l d2003bl) . 
i.e. using /3 = 5 (see Sect. 12.2.2b . 

3.2.1. Efficiency of the excitation 

A common feature of every modelling of stochastic excitation by 
turbulent convection is that the energy injecti on rate is locally 
proportional to the kine t ic ene rgy flux (e.g., iGoldreich et all 
1 1994 ISamadi & Goupill 1200 lb . The contribution due to the 
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Fig. 4. Top: Cumulative work integral integrals (W) versus the 
logarithm of the temperature for the n = 7, n = 9 and n 
12./ modes. An increases outward of W corresponds to 
a driving while a decrease inward, to a damping. The surface 
values give the final values of the work integral. By convention, 
the mode is found stable if the surface value of W is negative 
and unstable if positive. Bottom: Real part of the perturbation 
of pressure (6P/P) related to the n = 7, n = 9 and n = 12, 1 = 
modes versus the logarithm of the temperature. The horizontal 
dots lines permits to visualise the nodes of the eigen modes, and 
the vertical lines delimit the damping region. 



Fig. 5. Top: Rates (P) at which energy is supplied to the modes, 
computed as detailed in Sect. 12.2.21 for a 10M o main-sequence 
model (see Sect. 12. 11 1. Note that only the contribution of the iron 
convective region is included in P. Bottom: Square of the radial 
derivative of the radial component of the eigenfunction for radial 
orders n — 4, n — 6 and n = 8 that correspond in the top panel 
of the first maximum, minimum and the second maximum of the 
energy injection rates. The vertical lines delimit the location of 
the convective region induced by the iron opacity bump. 



Reynolds stresses can be simplified as dSamadill2009l) 

Pac [vj TkinA4 (10) 

where A is a characteristic length, c s is the sound speed, and Tvm 
is the specific kinetic energy flux defined as 

^kin=P« r 3 ms (11) 

with p the density, and M lms the root mean square velocity given 
by the MLT. 

The helium convective region is very inefficient in exciting 
modes since the kinetic energy flux is small (see Sect. 12.1.21 ). 
The more favourable convective regions, in this star, are the con- 
vective core and the iron opacity bump induced one. Excitation 
of radial acoustic modes by the iron convective region is found 
to be efficient for two reasons. First, the iron opacity bump is lo- 
cated deeper in the star compared to the helium bump. Both, den- 



sity and turbulent velocities are higher. Thus, the energy avail- 
able for p modes is higher. Second, the efficiency of the exci- 
tation depends on the involved time-scales, i.e. the convective 
time-scale and the modal period. The latter is several hours and, 
using mixing-length arguments (i.e. r ~ A/u m \ t ), we find the 
convective time-scale to be also of several hours. Hence, excita- 
tion is nearly resonant. 

3.2.2. Periodic variation and seismic diagnostic 

The power supplied to the mode by the iron convection zone 
exhibits periodic variations with frequency (Fig. [5] top). Such 
a behaviour is explained by the location of the iron convection 
zone compared to the radial component of the displacement. As 
shown by Eq. (O, the energy injection rate depends on the radial 
derivative of the radial component of the displacement in the 
driving regions. Hence, the location of the displacement nodes 
relatively to the iron convective zone determines its derivative, 
subsequently the efficiency of the excitation. Such an oscilla- 
tion is therefore a signature of the excitation by the iron con- 
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vective region. Observational evidence would lead to associating 
the driving with this region since, as will be shown in Sect. |4] no 
similar results are obtained if we assume that the convective core 
drives the modes. 

The location of the iron convective zone can be inferred from 
the energy injection rates by identifying the radial order of the 
first maximum, which corresponds to the mode whose last node 
is located at the bottom of the convective region (see Fig. [5] bot- 
tom). A second way to proceed is to use the frequencies of two 
the maxima. The successive maxima of the energy injection rates 
satisfy the relation 



r 

-)r h 



k r dr' + cf> = mn 



(12) 



where m is an integer, k r the radial wave-number, (p the phase, R 
is the star radius, r\, is the location of the node which we iden- 
tify as an upper limit of the depth of the iron convective region. 
Hence, by considering the frequencies of two maxima of the en- 
ergy injection rates (denoted by the subscripts 1 and 2), the loca- 
tion of the bottom of the iron convective region can be deduced 
from the integral expression 



r-R r-R 

I k r2 dr' - I k rl dr' 

J n, J rh 



(13) 



From Eq. (I131 l. one immediately sees that only the frequen- 
cies of two consecutive maxima are needed to determine the 
depth of the bottom of the iron convective region. To illustrate 
this, we consider the limit of pure radial acoustic waves, such 
that Eq. ( fT3] l becomes 



(a>2 -Ml) 



Jn, c s 



(14) 



Furthermore, in the very crude simplification where the sound 
speed is constant, one obtains the very simple analytic relation 



R-r b 



{u>\ - u> 2 ) 



(15) 



Using the benchmark model described in Sect. 12.11 one has 
c.5 ~ 40 km.s" 1 and from Fig. [5] we have co\ as 150ji/Hz and 
u 2 ~ 250 a/Hz. From Eq. Q3]) we get (R - r b )/R x 0.044, while 
from Fig.[3]one has {R-rb) I R = 0.04. This order of estimate then 
illustrates that the detection of solar-like oscillations excited by 
the iron convective region could give information on the struc- 
ture of the iron convective region. 

To illustrate this point, figure [6] (top) displays the energy in- 
jection rates as function of the radial order for two models at dif- 
ferent evolution stages. The more evolved the star is, the deeper 
the iron convective zone is located since the effective tempera- 
ture decreases as the star evolves. It turns out that the modes for 
which the energy injection rate is maximum shift toward lower 
radial-order. As a results, it is possible to follow and localise the 
depth of the convective region (see Fig.|6j bottom). Note that en- 
ergy injection rates are sensitive to the location of the bottom of 
the convective zone to less than one percent of the star radius, as 
shown in Fig. [6] 

3.3. Mode amplitude 

Figure [7] displays mode amplitudes excited by the iron convec- 
tive region. It appears that using the same values of the parame- 
ters as in the Sun, the amplitudes are found above the solar max- 
imum and the CoRoT detection threshold. In addition, Fig. [7] 




5 10 
radial order 




3 4 5 6 7 8 

radial order of the first maximum of P 

Fig. 6. Top: Normalised energy injection rates for two evolu- 
tionary stages (the two first successive models in Fig. [TJ as a 
function of the mode radial order. Bottom: Location of the bot- 
tom of the iron convection zone, with respect to the surface, ver- 
sus the radial-order of the first maximum of the energy injection 
rates. Each point corresponds to the 10 M Q modes described in 
Sect. 12. H at different evolutionary stages displayed in Fig.Q] 



also shows that those results are very sensitive to the convective 
velocities. It shows that a factor two in the velocities leads to 
important discrepancies for mode amplitudes, so does the mix- 
ing length, since the convective velocities depend on the mixing- 
lentgh to the third power for inefficient convection. 

The development of numerical simulations is then the only 
way to obtain more reliable values for mode amplitudes. 



4. Driving by core convection 

In contrast with the uppermost convective layers of the star, the 
inner convective core transports a non-negligible part of the total 
energy flux; convection is highly efficient. Hence, one cannot in- 
fer properties of those regions by simply looking for similarities 
with the solar convective region. Indeed, the physics is differ- 
ent and the dynamical properties, which are of interest to com- 
pute mode amplitude, are very uncertain. For instance, in the 
adiabatic regime the mixing-length has no significant impact on 
the determination of the star convective flux. The typical length- 
scale and the convective velocities are, however, of crucial im- 
portance to determine the efficiency of the driving. 
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Fig. 7. Mode amplitude as a function of mode frequency, com- 
puted as described in Sect. I2.2.2l and l2~2l Diamonds correspond 
to the computation of the energy injection rates using /3 = 5 
and the velocities given by the MLT, while for stars the veloc- 
ities have been multiplied by a factor two and divided by two 
for squares (note that a factor two corresponds to a variation of 
30% of the mixing-length). The solid horiz ontal line correspond s 
to the maximum for the Sun as derived bv lMichel et al.l (I2009I) . 
and the dashed horizontal line to the CoRoT detection threshold. 




1 50 200 

v (/LiHz) 



300 



Fig. 8. Rates CP) at which energy is supplied to the modes, com- 
puted as detailed in Sect. I2.2.21 for our benchmark model, by the 
convective core. 



Recent 3D numerical simulations can give some hints on 
those properties for massive stars. A recent st udy of convection 
in the efficient regime has been proposed by iMeakin & Arnettl 
(I2007h using compressible 3D numerical simulations. They in- 
vestigated convective regions during oxygen shell burning and 
during hydrogen core burning for a 23 M star and proposed an 
extended description of the convective properties and compar- 
ison with mixing-length theory as well as other 3D numerical 
simulations of convective outer layers. It turns out that both the 
velocities and the typical-length-scale are in rather good agree- 
ment with mixing-length-theory with a parameter (a) of around 
1.7. 



4. 1 . Energy injection rates 

As already pointed out in Sect. I2.1.21 the kinetic energy flux in 
the convective core is an order of magnitude higher than in the 
iron convective region. Such a large flux is due to the large den- 
sities and then makes the convective core a good candidate to 
excite modes. However, the excitation efficiency depends also 
on the matching between convective time-scales and mode peri- 
ods and on the shape of eigenfunctions (see Eqs.[8]and|9]l. 

Figure. [8]displays the normalised energy injection rates as a 
function of frequency for radial p modes. One can distinguish 
between two regions, namely low v(e [70; 200] //Hz) and high 
v(e [200; 300] /iHz) frequencies. In the low frequency range, 
the decrease of the energy injection rates is dominated by the 
decrease of the ratio between the convective turn-over time scale 
and the modal period, and also the shape of the eigenfunctions 
since V depends on its radial derivative. In contrast, modes with 
higher frequencies exhibit nodes near the interface between the 
convective and radiative zones explaining the modulation of en- 
ergy injection rates. In terms of absolute values for the energy in- 
jection rates, as for the iron convection zone, they depend mainly 
on the injection length-scale but are not very sensitive to the 
mixing-length since convection is efficient (see Sect. I4.3I ). 

4.2. Diagnostic on the eddy-time correlation function 

Convective turn-over time-scales, evaluated in the MLT frame, 
are found to be of the order of a month, while the modal pe- 
riods of acoustic standing waves are several hours. Hence, the 
excitation will be in an off resonance regime and will therefore 
crucially depends on the way eddies are time-correlated. In this 
particular regime, the way the eddies are temporally-correlated 
can lead to order of magnitude differences in terms of mode am- 
plitude s, as shown bvlSamadi et al.ld2003al) for the solar p modes 
and bv lBelkacem et al.l d2009al) for solar g modes, where the ef- 
fect is much stronger. 

Therefore, we consider different modelling of this time- 
correlation function (see Xk m Eq. ([H)). Exponential and 
Lorentzian modelling reads 



Xk{u) = — e-l-^l 



noJk 1 + (u)/u> k r 
with the normalisation condition 



f 



Xk(oj)dcu = 1 



where cjk is the linewidth, defined as 



(Ok = 



kuk 



(16) 
(17) 

(18) 

(19) 



where A is a parame ter introduced to account for the uncertain- 
ties in defining qj k dMmforthj Il992t ISamadi & GoupilL 120011: 
IChaplin et al.Ll200"5h . the velocity Uk of the eddy with wav enum- 
ber k is related to the kinetic energy spectrum E(k) by (ISteinl 
119671) 

-2k 



r 



AkE{k) 



(20) 



The choice of an exponential function is motivated by exper- 
imental studies of homogeneous and isotropic turbulence. For 
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Fig. 9. Top: Analytical eddy-time correlation function (xk) as 
a function of the ration a>/a>k, where a>k is the convective fre- 
quency at the wave-number k. The vertical line corresponds to 
the ratio <o/o)k for the most energetic eddies and for the mode n = 
2. Hence for higher frequencies, the ratio u>/u>k will increase, 
and the driving will be in the off resonance regime. Bottom:. 
Normalized mode amplitude computed using the eddy-time cor- 
relation functions as described in the upper panel. 



haviour of the mode amplitudes with frequency. The shape of the 
eddy-time correlation function has a great impact on the shape 
of mode amplitudes. A Lorentzian function, which slowly de- 
creases with frequency, results in a relatively small decrease of 
mode amplitude with increasing frequency. In contrast, an expo- 
nential Xk leads to a more steep slope for mode amplitude that 
makes the detection less likely. 

In conclusion, we have shown that the efficiency of the driv- 
ing by core convection depends crucially on the way the eddies 
are temporally-correlated. It is worthwhile to note that only a 
Lorentzian^ leads to an efficient driving and that the resulting 
mode amplitudes decrease with frequency in contrast with the 
driving by the iron convection zone. Hence, it constitutes a seis- 
mic diagnostic to identify the driving by the convective core and 
also to determine the way the eddies are temporally-correlated, 
an essential constraint on the dynamical properties of turbulent 
convection in cores of massive stars. 



4.3. Mode amplitudes 

Mode amplitudes are presented in figure [TOj To investigate the 
sensitivity to the injection length-scale, we compute amplitude 
with an upper limit that is the convective core size and a lower 
limit, the minimum of the pressure length-scale in the convective 
core. In the more optimistic estimate, relative magnitudes are 
found to reach up to 90 ppm for the lowest frequency mode and 
tens of ppm for higher frequency modes. Such a result demon- 
strates that excitation by the convective core can lead to am- 
plitudes well above the CoRoT detection threshold 1 ppm) 
as well as the solar maximum 3 ppm). In the pessimistic 
case, the amplitude is near the CoRoT threshold, which makes 
the detection more difficult. Note also, that those results assume 
a Lorentzian description of the eddy-time correlation function 
(Xk), and that using another prescription such as an exponential 
decrease or a Gaussian one lead to very small amplitudes well 
below one ppm. 

Nevertheless, it is worthwhile to note that the frequency be- 
haviour associated with the excitation by the convective core is 
very different from that obtained for the excitation by the iron 
convection zone. Hence, it constitutes a seismic signature that 
would be a powerful tool to identify the driving region. 



instance. iMordant et al.l d2004l) have demonstrated that the time- 
correlation of lagrangian velocities follows an exponential de- 
crease in agreement with the Kolmogorov (1941) phenomenol- 
ogy. However, the solar turbulent convection exhibits different 
physical conditions compared to experimental studies. The very 
large Reynolds numbers as well as the presence of coherent 
large-scale structures (plumes) make it likely that a specific de- 
scription of Xk is needed. Indeed, a Lorentzian description can 
be ad opted and has proven to better reproduce both the observa- 
tions (Belkace m et all [2006b|) and the 3D numerical simulations 
for bo th p and g modes dSamadi et"ail l2003ah iBelkacem et al.L 
I2009al) . 

Figure [9] (top) displays analytical^ versus the ratio u>/u>k, 
where o>k is the eddies characteristic frequency at the wave- 
number k. It turns out that by considering the most energetic 
eddies, excitation of p modes will occur in the off resonance 
regime with a ratio co/iOk higher than around 20. In that regime, 
it appears that the choice of Xk is essential as well as the width 
Uk- Note also that a Gaussian modelling leads to a vanishing 
driving since for such ratios u>/(L>k, the correlation is almost zero. 
Figure [9] (bottom) also shows the influence of the frequency be- 



5. Concluding remarks 

We performed an exploratory study of the driving (and damp- 
ing) of acoustic modes by turbulent convection in a 1OM star. 
It turns out that both the convective region associated with the 
iron opacity bump and the convective core are able to efficiently 
drive acoustic modes so that mode amplitudes are found above 
the CoRoT detection threshold 1 ppm) as well as the solar 
maximum (^ 3 ppm). 

However, uncertainties associated with the computation of 
mode amplitudes are important and are related to the poor 
knowledge of the dynamical properties of those convective re- 
gions. For the iron convective region, the main uncertainties 
come from the convective velocities as well as the injection- 
lentgh scale. For the convective core, the way eddies are 
temporally-correlated is crucial since excitation occurs in the off 
resonance regime. 

Computation of mode amplitudes also permitted us to em- 
phasise potential seismic diagnostics. First, driving by the iron 
convective region results in an oscillation of the energy injec- 
tion rates (and mode amplitudes) versus mode frequencies. It 
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Fig. 10. Mode amplitude versus frequency, computed using the 
benchmark model (see Sect. 12.11 ) and a Lorentzian function to 
describe the eddy-time correlation function (see Eq. (TTTIi). The 
upper line (with diamond dots) corresponds to the upper limit 
for the injection length-scale that is. The lower line (with stars) 
corresponds to an injection length-scale that corresponds to the 
minimum of the pressure scale height in the convective core. 



has been shown that this behaviour is potentially useful to deter- 
mine the extent and the depth of this convective region, as well 
as the dynamical properties of the driving (e.g., turbulent ve- 
locities, injection length-scales). Second, the way the eddies are 
temporally-correlated will determine the efficiency of the driv- 
ing by the convective core and the frequency behaviour. Hence, 
it constitutes an opportunity to strongly constrain^. 

The detecti on of solar-like oscillati ons in a massive star such 
as VI 449 Aql dBelkacem et aUl2009bb and the determination of 
the mode parameter in this star are a promising way to iden- 
tify the driving region, since both have very different signa- 
tures, and would give physical constraints to the still poorly 
known treatment of convecti ve regions. The prese nce of high 
frequency power in SPB stars (iDegroote et al.l 120091) would also 
help to constrain the nature of the driving zone since these stars 
have luminosities definitively below the threshold value for the 
presence of a convec tive zone associated with the iron peak 
dCantiello et al.L l2009t) . It should however be kept in mind that 
microscopic diffusion coupled with radiation acceleration could 
induce an accumulation of iron in the iron peak. A convective 
iron zone could thus also be present in stars whose luminosity 
are slightly below the threshold. 

Acknowledgements. We are indebted to J. Leibacher for his careful reading of 
the manuscript and to R. Samadi and M.J. Goupil for their helpful remarks. K.B. 
also acknowledges financial support through a postdoctoral fellowship from the 
"Subside federal pour la recherche 2009", University of Liege. 



Cantiello, M., Langer, N., Brott, I., et al. 2009, ArXiv e-prints 
Chaplin, W. J., Houdek, G., Elsworth, Y., et al. 2005, MNRAS, 360, 859 
Christensen-Dalsgaard, J. & Dappen, W. 1992, A&A Rev., 4, 267 
de Ridder, J., Barban, C, Baudin, R, et al. 2009, Nature, 459, 398 
Degroote, P., Aerts, C, Ollivier, M., et al. 2009, ArXiv e-prints 
Dupret, M.-A. 2002, Bull. Soc. Roy. Sc. Liege, 5-6, 249 

Dupret, M.-A., De Ridder, J., Neuforge, C., Aerts, C., & Scuflaire, R. 2002, 
A&A, 385, 563 

Dupret, M.-A., Thoul, A., Scuflaire, R., et al. 2004, A&A, 415, 251 
Dziemblowski, W. 1977, Acta Astronomica, 27, 95 

Dziembowski, W. A. 2007, Communications in Asteroseismology, 150, 175 
Dziembowski, W. A., Moskalik, P., & Pamyatnykh, A. A. 1993, MNRAS, 265, 
588 

Dziembowski, W. A. & Pamiatnykh, A. A. 1993, MNRAS, 262, 204 

Dziembowski, W. A. & Pamyatnykh, A. A. 2008, MNRAS, 385, 2061 

Gautschy, A. & Saio, H. 1995, ARA&A, 33, 75 

Gautschy, A. & Saio, H. 1996, ARA&A, 34, 551 

Goldreich, P., Murray, N., & Kumar, P. 1994, ApJ, 424, 466 

Grevesse, N. & Noels, A. 1993, in Origin and Evolution of the Elements, ed. 

N. Prantzos, E. Vangioni-Flam, & M. Casse, 15-25 
Grigahcene, A., Dupret, M.-A., Gabriel, M., Garrido, R., & Scuflaire, R. 2005, 

A&A, 434, 1055 
Iglesias, C. A. & Rogers, F. J. 1996, ApJ, 464, 943 
Iglesias, C. A., Rogers, F. J., & Wilson, B. G. 1992, ApJ, 397, 717 
Kippenhahn, R. & Weigert, A. 1990, Stellar Structure and Evolution (Springer- 

Verlag (Berlin)) 

Lighthill, M. J. 1952, Royal Society of London Proceedings Series A, 211, 564 

Meakin, C. A. & Arnett, D. 2007, ApJ, 667, 448 

Michel, E., Baglin, A., Auvergne, M., et al. 2008, Science, 322, 558 

Michel, E., Samadi, R., Baudin, F., et al. 2009, A&A, 495, 979 

Mordant, N., Leveque, E., & Pinton, J.-F. 2004, New Journal of Physics, 6, 1 16 

Moskalik, P. & Dziembowski, W. A. 1992, A&A, 256, L5 

Pamyatnykh, A. A. 1999, Acta Astronomica, 49, 119 

Pamyatnykh, A. A., Handler, G., & Dziembowski, W. A. 2004, MNRAS, 350, 
1022 

Pesnell, W. D. 1990, ApJ, 363, 227 

Samadi, R. 2009, Lectures Notes in Physics, In press 

Samadi, R., Belkacem, K., Goupil, M. J., Dupret, M.-A., & Kupka, F. 2008, 

A&A, 489, 291 
Samadi, R. & Goupil, M. . 2001, A&A, 370, 136 

Samadi, R., Nordlund, A., Stein, R. R, Goupil, M. J., & Roxburgh, I. 2003a, 
A&A, 404, 1129 

Samadi, R., Nordlund, A., Stein, R. F, Goupil, M. J., & Roxburgh, I. 2003b, 
A&A, 403, 303 

Scuflaire, R., Theado, S., Montalban, J., et al. 2008, Ap&SS, 316, 83 
Stein, R. F. 1967, Solar Physics, 2, 385 

Thoul, A., Aerts, C, Dupret, M. A., et al. 2003, A&A, 406, 287 



References 

Aerts, C, Thoul, A., Daszynska, J., et al. 2003, Science, 300, 1926 
Alexander, D. R. & Ferguson, J. W. 1994, ApJ, 437, 879 
Balmforth, N. J. 1992, MNRAS, 255, 603 

Bedding, T. R. & Kjeldsen, H. 2007, Communications in Asteroseismology, 150, 
106 

Belkacem, K., Samadi, R., Goupil, M.-J., & Dupret, M.-A. 2008, A&A, 478, 163 
Belkacem, K., Samadi, R., Goupil, M. J., et al. 2009a, A&A, 494, 191 
Belkacem, K., Samadi, R., Goupil, M. J., & Kupka, F. 2006a, A&A, 460, 173 
Belkacem, K, Samadi, R., Goupil, M. J., Kupka, R, & Baudin, F. 2006b, A&A, 
460, 183 

Belkacem, K, Samadi, R., Goupil, M.-J., et al. 2009b, Science, 324, 1540 



